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Abstract. The two Fresnel Integrals are real and imaginary part of the in- 
tegral over complex-valued exp(ix 2 ) as a function of the upper limit. They 
are special cases of the integrals over x m exp(ix n ) for integer powers m and n, 
which are essentially Incomplete Gamma Functions. We generalize one step 
further and focus on evaluation of the integrals with kernel p(x) exp[i(j>(x)] 
and polynomials p and <j>. Series reversion of <j> seems not to help much, but 
repeated partial integration leads to a first order differential equation for an 
auxiliary oscillating function which allows to fuse the integrals and their com- 
plementary integrals. 



1. Introduction 

1.1. Motivation. A criticism to the orthogonal basis of the Zernike Polynomials 
on the unit circle is the lack of the minimax property-if compared for example to 
the Chebyshev Polynomials on the one-dimensional interval [25]. The radial poly- 
nomials of higher order, starting already with Zf 1 , Z\ and Z^ 2 , attain amplitudes 
at the local maxima and minima within the circle which are smaller (in absolute 
value) than the amplitude at the rim. One way to equilibrate these, to enhance the 
extrema in the interior and to weaken the extrema along the perimeter, is to add 
higher powers of the radial distance to the polynomials. An alternative, which trig- 
gered this work, may employ sines and cosines of polynomials of the radial distance 
which evidently oscillate up and down with constant amplitude. The requirement 
of orthogonality and normalization then puts focus on radial integrals of the form 

(1.1) W«) = / p(x)e l ^ x) d Xl 

Jo 

where <f> is the sum or difference of two phase-polynomials of the basis. 

The concept applies to Zernike Polynomials in higher dimensions as well, so p are 
(small) integer powers of the radial coordinate, representing the radial part of the 
Jacobian determinant of the transformation from Cartesian to Hypcr-sphcrical Co- 
ordinates. (We assume that the angular variables are represented by Hyperspherical 
Harmonics, which ensures that these oscillations oscillate with amplitudes of equal 
magnitude already by design.) This work deals with the evaluation of the integrals 
I p ,4> given a finite upper limit u plus the polynomials p and cj>. The adjustment 
(variation) of the coefficients of a set of polynomials <j> to achieve orthogonality and 
to establish an explicit basis of this type is not considered here. 
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1.2. Overview. Section 1.3 summarizes some elementary closed-form solutions. 
Section 2 is an overview of known results concerning the cases with constant p{x) 
and <j)(x) being a power of x, which includes the Fresnel Integrals eo ipso. Section 
3 widens the perspective to the cases where both p(x) and <p(x) are integer powers 
of x. The evaluation is covered by the standard theory of the Incomplete Gamma 
Function. Section 4 looks at series reversions if 4>{x) are polynomials, leading to 
barren results. Section 5 shows that asymptotic expression for the complemen- 
tary integral with one limit bound at infinity is obtained easily via an auxiliary 
inhomogeneous differential equation. The fusion of the series representations of 
that auxiliary function at small upper limits on one hand and its asymptotic series 
on the other is prosecuted by finding the adaptive free parameter of the general 
solution of the differential equation. 

1.3. Simple Cases. The simplest cases are polynomials <f> of zero or first degree 
[29, 2.635]: 



(1.2) / xcosaodx 

(1.3) J xsmaodx 
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(1.5) / x sin(o!o + a\x)dx = 
Repeated partial integration yields 

/ iax J™ „l 

with real and imaginary parts [29, 2.323,2.633] 

(1.7) / x n cos(ax)dx = > ^—r . sin(az + kn/2), 

J (n - k)\ a k+1 

J k=0 y ' 

/■sr-\ n\ x n ~ k 
x n sin(ax)dx = — > — — , , , cos(a.T + kir/2). 
y ' In- k)\ a k+1 
k=0 y ' 

(Erratum: A factor (— l) k is missing in the last term of [29, 2.633.6]). Further 
obvious special cases are 

(1.9) [ x^ 1 sin(x n )dx = -- cos(x n ); [ x 11 ' 1 cos(x n )dx = - sin(x n ); 

./ n J n 



(1.10) 



c™ 1 e lx "dx = — ( 
in 



Remark 1. The polynomial in the argument of the trigonometric function may be 
assumed to have minimum low degree 1, because its constant term is easily moved 
to the front of the integral with [1, 4.3.16,4.3.17] 

I I I 

(1.11) sin(^^ ctjx-i) = sin qp cos(y~^ ajX-*) + cos otp sin(^~] ctjX-'), 

j=0 3=1 j=l 
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Figure 1. The integral J™ sm(x n )dx for three different powers n. 
Extrema are at u n = lir with I = 0, 1, 2, 



(1.12) 



i (y ] ajX-i) = cosao cos(^^ ctjX J ) — sinao sin(^^ ajX J ^ 

3 = 1 



3=0 



3 = 1 



which are real and imaginary part of 

I I 
(1.13) exp(i OijX-*) = exp(iao) cxp(i^^ otjX^). 



3=0 



3 = 1 



Combinations of an odd polynomial p{x) with an even polynomial <p(x) are re- 
duced in a first step by the substitution y = x 2 , 



(1.14) 



f u 1 /"" 2 

/ x 2m+1 exp{tx 2n )dx = - / y m exp(iy n 
Jo 2 Jo 



)dy. 



2. Fresnel Integrals 



The integrals J™ s'm(x 2 )dx and J™ cos(x 2 )dx define — up to some unimportant 
constant scale factor— the Fresnel Integrals [1, 7.3] [29, 8.25] [6, 30, 48, 8, 46]. In 
this section we look at the generalized integrals J sin(a; ra )<fa: and J cos(x n )dx for 
integer n = 2, 3, 4, 
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Figure 2. The integral J™ cos(x n )dx for three different powers 
Extrema are at u n = (I + 1/2)tt with I = 0,1,2, 



The substitution x = y 1 ^ 71 followed by iterated partial integration generates 
1 



(2.1) J e lx "dx = - J y 1 / n - l e tv dy 

= l y i/n-i e iy _ 1(1 _ n f y^-2 e *y dy 
in in n J 

= l y Vn-i e *_l(I_i)_ J_(I_i) tf V»-V»+ 4~ (--l)(--2) / y l/n ^e^dy 



The product of the factors (— — 1)(— — 2)(— — 3) • • • has factorial growth which 
induces asymptotic convergence for this representation. Because this is the well- 
known series of the confluent hypergeometric function in (3.2), we need not dwell 
into this further [1, 13.5.1]. 
For small x, the Taylor series 



(2.2) / muu-W,- = / J2 |L- — H'dx = £ 
J J i=o { 1=0 



r 2nl-\-n-\-l 



(2nl + n + l) (2Z + 1)! 



and companion series for cos(a;™) converge quickly and are the standard target 
for Chebyshev approximations in that area of numerical representation of special 
functions. 
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Figure 3. The integral J™ x 2 sm(x n )dx for three different values 
of n. 

Algorithm 1. [21 , p 63] If the coefficients b define a power series 



(2.3) 



s>0 



then the coefficients of the Neumann expansion 
(2.4) z"f(z)=J2asJv+s(z) 



s>0 



L-/2J 



(2.5) 



= (!/ + *) £ 2" 



+S -2mI> + S-m) 



6,-s 



Application of this algorithm with v = 1 + 1/n, z = x™ and n = 2 to Eq. 
yields the well established [29, 8.515] [5] 



(2.6) 



/ sm(x 2 )dx = J2s+z/2{x 2 ) 
J z s >0 



For n — 3 to 5, the Neumann series 

(2.7) fsm(x n )dx = d n J2tn, s J2s+„(x n ) 



s>0 




Figure 4. The integral J™ x 2 cos(x n )dx for three different values 
of n. 

Table 1. Table of Neumann series coefficients (2.7) for n = 3 with 
d 3 = 2 1 / 3 7rV3/r(2/3) w 5.062876576879227 



4/27 « .14814814814814814814814815 

1 136/729 « .18655692729766803840877915 

2 1120/6561 « .17070568510897729004724889 

3 31912/177147 w .18014417404754243650753329 

4 2491372/14348907 « .17362799828586247022159946 

5 23052640/129140163 .17850867975131795365629204 

6 608995904/3486784401 .17465831951793224739736353 

7 16739336224/94143178827 .17780721272181224941292368 

8 148413541000/847288609443 w .17516291302153317339294102 

9 12177019210000/68630377364883 .17742899977453386676382033 

10 325122763152640/1853020188851841 w .17545559681899144476569162 



are detailed in Tables 1-3 up to order s = 10. 

Remark 2. According to (3.8), lim s _>. 00 d n £,n, s = T( ) s i n (^")- ^ e observation 
that £„ jS siays practically constant as s grows tells us that this expansion converges 
poorly for medium and large upper limits of the integral. 
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Table 2. Tabic of Neumann series coefficients (2.7) for n = 4 with 
d 4 = 2 3 / 4 7r/T(3/4) ps 4.31 f 60f 09908f 8558898775f 



1/8 ps .12500000000000000000000000 
11/64 ps .17187500000000000000000000 
159/1024 ps .15527343750000000000000000 
1347/8192 ps .16442871093750000000000000 
41531/262144 ps .15842819213867187500000000 
341309/2097152 ps .16274881362915039062500000 
5350155/33554432 ps .15944704413414001464843750 
43506995/268435456 ps .16207618638873100280761719 
2747379155/17179869184 ps .15991851425496861338615417 
22228087705/137438953472 ps .16173062398593174293637276 
352241109985/2199023255552 ps .16018071163898639497347176 



Table 3. Tabic of Neumann series coefficients (2.7) for n = 5 with 
d 5 = 2 1 / 5 7r/[sin(7r/5)r(4/5)] ps 5.27349462002150700879064 



e 



5..S 



2/25 ps .080000000000000000000000000 
72/625 ps .11520000000000000000000000 
8118/78125 ps .10391040000000000000000000 
214544/1953125 ps .10984652800000000000000000 
5179592/48828125 ps .10607804416000000000000000 
663616816/6103515625 ps .10872697913344000000000000 
16287458776/152587890625 ps .10674148983439360000000000 
2065603429328/19073486328125 ps .10829710907555184640000000 
51039744111764/476837158203125 ps .10703810144347409612800000 
1288445436120032/11920928955078125 ps .10808263692999989395456000 
159738815063405788/1490116119384765625 ps .10719890415611250403704832 



3. Product with Powers, Incomplete Gamma Functions 

Introducing a power of the free variable in the integrand leads from (2.2) to the 
more general 

(3.1) 

/f 00 m+2nl+n 00 m+2nl+n+l 

m** - /E^ri-i* - g (m+ 2l7;, i+1 )-(2 l+1)! 

_,m+n+l / 1 1 m+1 ™2n\ 
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This is the imaginary part of the following confluent hypergeometric function [37] : 

/r 00 jl T m+nl 
x m exp{ix n )dx = J\ dx 



00 jl x m+nl 00 jl x m+nl+l 

^ (to + nl + 1) II 

1=0 1=0 y ' 



ro+1 

n 

1 1 " 1 V 1 + m + 1 



1*1 



TO 4 

These underivatives are illustrated in Figure 3 and 4. 

Remark 3. By the rules (a) 2r = (a) r (±±^) r 4 r and (a) 2r +i = a(±±^) r (l + § )A r for 
the Pochhammer Symbol [40] , the even and odd parts of a generalized hypergeometric 
function r F s (z) are [22] 
(3.3) 

ai Q2 l+a2 a r l+a r 

2' 2 '2' 2 ''''2' 2 



E (ai)n(Q2)n ■ • • (a r )n Z n _ 
IhA (h„\ ...(h\ ~ 2rI ' 2s + 1 \ hi L+El l+fe . . . br. 1+br I 



n=0,2,4,6,. 

and 



(3-4) £ 



(bi) n (b 2 )n ■ ■ ■ (b s )„ n\ 

(ai) n (a 2 ) n ■ ■ ■ (a r ) n z T 



2' 2 ' 2 ' 2' 2' 2 ' 2 



(bl)n(b 2 )n ■ ■ ■ (b s )n Tl\ 

n r a . f 1+ai 1 , 01 l+a 2 1 I _2_ ... 1 1 

5 7 _ Z 2 r^2s+l I -, fc^ l+b 2 , 63 1+br 1 i &r. 3 

j=l u j \ 2' ~ r 2'2' ~ r 2' 2 ' 2 ' 2 

where parentheses with subindex denote Pochhammer 's symbol, 

(3.5) (a)„ a(a + l)(a + 2) ■ ■ ■ (a + n - 1) = T(a + n)/T{a). 

Extraction of real and imaginary parts of hypergeometric functions with purely imag- 
inary argument z is an application of these formulas. 

Remark 4. The increase in amplitude observed in these plots for large x if 4> has 
low degree n does not hamper the numerical treatment. The partial integration 

m—n-\-l ix n 



(3.6) in j x m e lx dx = x m - n+1 e lx -(m-n+1 

decrements the exponent to in the integrand — closely related to a contiguous relation 
of \F\ in (3.2) [1, 13.4.7]. Combined with (1.10), the task that remains is to handle 
the cases to < n. More generally speaking, if the degree of the polynomial p{x) 
in the integrals of the form J p(x)e l ^ x ^dx is at least the degree of the polynomial 
4>(x) minus one, a polynomial division of p through the first derivative of 4>, p(x) = 
p(x)cj)'(x) + r(x), followed by partial integration of the term J p(x)(j)' (x)e l ^ x ^ dx 
reduces the task to the format J p' (x)e l ^ x ^ dx and J r{x)e l ^ x ^dx where polynomial 
degrees are degp' < deg^' and degr < deg^' — this reduction executed recursively 
if needed. 

The leading term in a Laurent expansion is [1, 13.5.1] 

(3.7) 1 F 1 ( 1 7t; +1 | ixA -+ !Ll±l T (Jl±ly^ +mx ^ m+ i)^ 

which is inserted into (3.2), 

(3.8) x m exp(ix n )dx = Ir(-_±-.)e-( m+1 )/( 2 "), 
./n n n 
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and by the scaling substitution ic 1 '" — > y, 
(3.9) 

x m cx P (icx n )dx = n ! w T<?^±y^ m+V >iw = i(!)( 1 +-)/»r(^^). 



j lc v_j.-rT»)/« n n c 

Remark 5. Particular imaginary parts of these limits are 
(3.10) 

f 0.626657068657750125603941, n = 2; 

/ sin(a;")da; = -T(l/n) sin ^- w i 0.446489755784624605609282, n = 3; 

[ 0.346865211023809496042035, n = 4. 

Remark 6. vl simple pattern emerges for the k-fold repeated integrals [20], 
(3.11) 

p p 00 

/ • • • / x m eMix n )dx = 



(m + nl + 1) (m + nl + 2) • • ■ (m + nl + k) l\ 

x m+k / m+1 m+2 _ _ _ ro+fc 

fc fc I 1 _i_ ra±l i _l »»+2 i , m+fc i ix 



(m + l)(m + 2) ■ ■ • (m + k) K K \ 1 + 1 + ^ 

Kummer's transformation of (3.2) is [1, 13.1.27] 

(3.12) J x m cxp(ix n )dx = V n , m (x)e l 

with 



m+l 



(3.13) v n , m {x) = —— 1 F 1 ( „ 

m + l \ 1 H — 

Real and imaginary part of (3.12) are 

(3.14) / x m cos(x n )dx = W n , m (x) cos(x n ) - QV n , m (x) sm(x n ); 



(3.15) / x m sm(x n )dx = UV n , m (x) sin(x") + SV n , m (x) cos{x n ). 

Remark 7. The same result could have been derived by repeated partial integration, 
the initial two steps being 

(3.16) [ x m e lx "dx = — !— -x m+1 e"™ — [ x m+n e ix "dx 

J m+l m+l J 

_ x ™+l e «" + _ _Z™ m+n+l e *x» + HlO! / 2 ,m+2n e «" d , 



m+l (m+l)(m + n+l) (m + l)(m + ra + 1) 

Gathering each second term of the hypergeometric series of (3.13) 

x m+l « (-)^(ffHl 



(3.17) 3tV n , m (x) 



^(l + ^W^ + l)! 



1 



(m + l)(l + ^)^ 2 V 1 + ^,1 + ^ ' 4 

A corresponding calculation of the other half of the terms yields 

1 / 1 r 2 " 

(3.18) s St V n , m (x) = —-x m+1 1 F 2 l ! , m+ A , m+1 I- — 
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-0.8 
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2 2.5 3 3.5 4 4.5 5 

x 



Figure 5. The value of iF 2 (l; \ + ^,1 + the main 

constituent of $lV n fl(x) in (3.18), for three different degrees n of 
the phase polynomial. 



Remark 8. One might hope that splitting off the factor exp(zx n ) in (3.12) removes 
the major oscillations and V n ^ m contains less wiggles. Figure 5 illustrates that this 
target is missed, and Section 5.3 puts this into a wider perspective. 

Remark 9. The continued fraction of (3.13) is [17, 37, 33] 



(3.19) ifW 1 



l+a 



_|_ 0!(2+q) | (l+a)(3+q) | l!(4+q) | 



(l+a) 



(l+q)(2+q) 
X 



(l+q)(2+q)(5+q) _ 2!(6+a) 

2! ' (l+a)(2+q)(3+q) 



(l+q)(2+q)(3+q)(7+q) . 3!(8+q) 

3! ~T (l+q)(2+q)(3+q)(4+q) 



(l+a)(2+q)(3+q)(4+q)(9+q) 
4! 



Fields and Wimp and Luke provide a series expansion in terms of Bessel func- 
tions [24, (2.14)] [36]. 

Eq. (3.2) transforms into an Incomplete Gamma Function [1, 13.6.10] [29, 2.632] [3, 
26, 27, 10, 2, 4]: 



(3.20) 



x m exp(ix n )dx = Ii( m+1 )/" 7 (^-i-,-^ n ), 
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which is sometimes phrased as [1, 6.5.3] 

(3.21) l(a,z) =T(o) -T(a,z). 

This association opens a wealth of literature to numerical evaluation. A partic- 
ularly nicely converging continued fraction representation is [17, 12.6.34] [42] 

(3.22) r M = ^L! t^l 2( S -2) 3(-3) 4(« - 4) 

1 + z - s+ 3 + z - s+ 5 + z - s+ 7 + z - s+ 9 + z - s+ 

Remark 10. This can apparently be re-binned in terms of associated Laguerre 
polynomials as [16] [43, (41)] 

(1 - s) n 1 



(3.23) T(s,x) = x s e' x Y^ 



4. Series Reversion 

The final chapters of the manuscript are concerned with phases 4>(x) which are 
not integer powers but full-fledged polynomials of x. 

4.1. Expansion Around the Origin. If the leading (lower) degree of the phase 
polynomial is 1 (so ceo = after use of (1-11) and (1-12) and a± ^ 0), a substitution 

i 

(4.1) y = Y, a o xJ 

3=1 

in the integrand is available, associated with a series reversion of the polynomial, 

oo 

(4.2) x = J2&y j > 

i=i 
i=i 

The idea is that preserving the sine/cosine/exponential in the integrand allows to 
deal efficiently with the oscillations in the spirit of Filon quadratures. 



(4.3) dx/dy = 



Algorithm 2. The intermediate task is to determine the (infinite) set of the /3j 
from the (finite set of) otj. The (3j are tabulated as a function of the ctj by 
Abramowitz and Stegun [1, 3.6.25] for j < 7, and by Orstrand [38] for j < 13. 
In the general case, equalize the coefficients of equal powers of y in (4-l)> which is 

l I oo \3 

(4.4) y = J> < X>V* 
j=i \fc=i / 

/ oo / ■ \ 

*=J tx + 2t 2 -\ h mt m = t 

t l +t 2 ^ h t m = j 

t\,t%,..., t m > 

involving the multinomial coefficients in parentheses [39, 38, 11, 23, 47, 19]. The 
maximum index m that contributes with t m > to the sums iU = t and 
^2™ ti = j is m = t — j + 1. The number of multinomial terms can be read off 
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entry A008284 in the Online Encyclopedia of Integer Sequences (OEIS) [41]. Their 
evaluation may be avoided by employing a convolutional recurrence [29, 0.314]. 
In the order 0(y 1 ), 

(4.5) l = c*i/3i. 

In the orders 0(y ), t > 1 ; fit is derived by isolating fit in the term with j = 1 in 

min(t,Z) , v 

( 4 - 6 ) °= E °y E ( h J... trr )##-#r 



j=1 *i + 2* 2 + • • • + mt m = t 

h + 1 2 + ■ ■ ■ + t m = j 

tl, t 2 , ■ ■ ■ , t m > 



£ (oo.!.oi 



ti + 2t- 2 + ■ ■ ■ + mt m = t 
*i + * 2 + • • • + t m = 1 
ti, t 2l ■ ■ ■ , t m > 

min(i.Z) . * 

+ E ** E UJ..O##"-^' t>L 

J - 2 ti + 2t 2 H h m* m = t 

ti + 1 2 H h * TO = j 

ti,t 2 , . . . ,t m > 

which yields the recursive algorithm 
(4.7) 

min(i,Z) , 



«iA = - E «i E ( M2 J .. tm ^-^>i. 

J - 2 *i + 2* 2 + ■ • • + m* m = * 
*1 + * 2 + ' ' ' + t m = j 
ti, *2j • • • j *m > 

In practise, the problem is solved by a perturbative expansion, starting from the es- 
timate a?*- 1 -* = y/ai and inserting iteratively previous estimates into the constituent 
equation: = y — J2]= 2 ctj(x f - s ~ 1 ^y [7]. 

For a set of simple polynomials, the [3 coefficients are gathered in Appendix A. 
The application of the series reversion linearizes the phase argument in the 
trigonometric function. 

Example 1. 

„ l „ oo oo 

(4.8) J xexpii^ajx^dx = J kf3 k y k ' 1 ){^2 PjV 5 ) cxp(iy)dy, 



j=l " fc=l 3=1 



oo . 

E 7t / y* ex p(^) d y' 



1=1 

where the series coefficients 7 are derived from convolution type products 

t 

(4.9) lt = J2Wt-k+v 

k=l 
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Approximation of the representation of the series by truncation at some maxi- 
mum index of the (3 reduces the problem for each term in (4.8) to the closed forms 
(1.8). 

Remark 11. // the polynomial in the phase argument has degree I = 2, the linear 
term in the polynomial «o + cti% + O-2X 2 may be eliminated by the substitution 
z = x + ai-i/(lai). The remaining absolute term can also be removed as mentioned 
in Remark 1, so J cos(ao + aix + a>2X 2 )dx and J sin(ao + ot\x + ot2X 2 )dx have closed 
form expressions in terms of Fresnel Integrals [29, 2.549.3, 2.549.4]. 

4.2. Local Series. The main reason not to use the approach of Section 4.1 in 
numerical application is that the series reversion is only unique/stable up to the 
argument of the first maximum of the phase polynomial y(x). The administrative 
problem is to monitor branch cuts of the inverse x(y) and to monitor the stability 
of the results as a function of the truncation order of the j3. The quick growth of 
the sequences in Appendix A demonstrates the problems of convergence, and for 
small upper limits within the radius of convergence, the competing Taylor series 
are simpler to handle. One improvement is to split the interval of integration and 
to revert Taylor series of the phase polynomial anchored in the middle of the sub- 
intervals. 

The power series expansion of the integrand near some x is found with the 
binomial expansion 



(4.10) exp(i(x + e) n ) = exp(ix n + iJ2 (?V 



1=1 



1-1 



cxp(ix n )exp \ [l) e 



1=1 



Assuming e is small, reversion of the polynomial in the argument on the right hand 
side generates the power series. We illustrate this procedure in the next three 
subsections. 

4.2.1. Squares. For quadratic chirp polynomials, n = 2, 

(4.11) exp ^i ^2 (^j <^ n ~^ = cxp(2zxe + it 2 ) = exp(2ixy), 
where 

(4.12) y^ e+ h 2 

is substituted in the integral. Series reversion 

(4 - 13) e = y -i y2 + (i^ y3 '(ir y4 + (^ y5 ---- 

with coefficients from the first row of Table 5 leads to the derivative 

(4,4) */* = 1 - ±y + ±±»> - ±±y° + ^ - . . . . g ^(-yf. 
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Integration over small e-environments near some x are then evaluated with 

(4.15) J exp(2ixe + ie 2 )de = £ J (_„)» exp(2ixy)dy 

— eventually with (1.8) — . In practise, these quadratic polynomials would rather 
be evaluated with the closed form procedure of Remark 11. 

4.2.2. Cubes. For cubic chirp polynomials, application of (4.10) reads 

(4.16) exp(i(x + e) 3 ) = cxp(ix 3 ) exp(3iex 2 + 3ie 2 x + ie 3 ) 
with one factor 

(4.17) exp(3ie.T 2 + 3ie 2 x + ie 3 ) = exp(3zx 2 (e + -e 2 + J-^e 3 )) = exp(3ix 2 y) 

x 3x z 

after substituting 

<« 8 > 

This reverts to [41, A218540] 
(4.19) 

_B_ 3 _10yl 4+ ^_ 5 _i5V3 J122_ 7 _280V3 
y (v^) 22 ' (VSx) 3 ^ (v^x)^ (v^z) 5 y + ^3xf V W3xy V 

E K i ,.2.7 + 1 _ \ " ^jV 3 2,7+2 
(V3a;) 2 J * ^ (VSa;) 2 ^ 1 U ' 

Remark 12. Alternatively, by canceling some square roots of 3, 

1 2 15 3 90 4 594 5 4158 „ 
4.20 e = y - -y 2 + j^y 3 ~ j^V + 7^4 V - +■•• 

x {3x) z {3x) d (3xp (3a;) 

with numerator coefficients tabulated in the OEIS sequence A025748 [41]. 

On the r.h.s., the coefficients ko = 1, K\ = 5, K2 = 66, K3 = 1122 etc. established 
in the numerators of odd powers of y obey n(2n + V)n n = 3(6?i — l)(3n — 2)re„_i. 
The recurrence is solved by 

(4.21) n n = 27 n (5/6)„ (l/3) n /[n!(3/2)„] 
with generating function 

(4.22) 5>«z"=2-Fi( 5/ 3/2 /3 |27z 

The coefficients Ao = 1, Ai = 10, A2 = 154, A3 = 2805,... in the numerators of 
even powers of y obey (n + l)(2n + 1)A„ = 3(6n — l)(3n + l)A„_i, such that 

,, 9 ^ . 27"(5/6)„ (4/3)„ l + 3n 

(n+l)!(3/2)„ l + n 

with generating function 

(4-24) ^ ^ = 3F 2 ( 5/ 3 6 ^A 1 1 27z 



n>0 
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The substitution actually calls for the derivatives 

(4.25) de/dy = V ^ ^ - V (2 / + 2 ^V 1 . 

j>0 vv ; j>0 vv ' 

So the cubic analog of (4.15) is 
(4.26) 

J eM^+^x^)de = J dyeM3™ 2 y)?£ ^Jte)*^ ^ 

for use in the right hand side of (4.10). 

4.2.3. Quartics. The quartic chirp with (4.10) reads 

(4.27) exp(i(x + e) 4 ) = exp(ix 4 ) exp(4iex 3 + 6ie 2 x 2 + Aie 3 x + ie 4 ) 

with small y in 
(4.28) 

Q 1 1 

exp(4ix 3 y) 



Aix 3 (e + ^-e 2 + -U 3 + -r^e 4 ) 

2x x 2 4x 3 



exp(Aiex 3 +6ie 2 x 2 +4ie 3 x+ie 4 ) = cxp 
substituting 

This reverts to 
(4.30) 

3 2 14 o 77 , 462 «- 2926 fi 19228 7 129789 8 

2 i x 2^ar 2 J a; cl 2 4 x 4 2 a a; 2°a; t) 2'a;' 



(2x) 

j>i v ' 



where the integer sequence rj is built up by [41, A048779] 



(4.31) JT ?j +2(5-4 J > i _ 1 =0; = ^— ^— ^- 



j>i 

With derivatives 

v (j + 1 )%+i 
1 {2x)i 



(4.32) de/dy - 



j>0 

the quartic variant of (4.15) is 
(4.33) 

/ exp(4iea; 3 + 6ie 2 x 2 + Aie 3 x + ie 4 )de = [ dyexp(4ix 3 y) V ^ % ^ Vi+1 (-yY. 

5. Associated Linear Differential Equation 

5.1. Perturbative. The iterated partial integration of Remark 7 defines a gener- 
alization of (3.12) to polynomials p(x) and <j>(x), 
(5.1) 

Ip.^x) = j p(x)e^dx = qi (x)e 1 ^ - I iq^'e^dx = e ^ ^ q 2k+1 (x), 

k—0 
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where 
(5.2) 

q (x)=p(x), q 1 (x) = I q (x)dx, q 2 k+i(x) = -i I q 2 k-i(x)(j>' (x)dx, k>0 



construct a sequence of polynomials of increasing degree. The sum Q2k+i equals 
V n ,m if P and 4> are integer powers of x; in that sense q(x) constructs the Kummcr 
transformation of I p ^(x). The sum over all <72fc+i is a perturbative construction of 
a solution q{x) to the first-order linear inhomogeneous differential equation 

(5.3) p(x)=iq(x) ( f>'(x)+q , (x), 
respectively integral equation 

(5.4) q(x) = / p(x)dx — / iq(x)(f>' (x)dx. 

Jo Jo 

It starts with the estimate that iq<f>' on the right hand side vanishes, which leads to 
the approximation that q is the integral over p. This is plugged into the first term 
of the right hand side, and solved by integration over p — iq<p' and so on. We are 
mainly interested in singular solutions with lower limit zero, so we insert lower 
limits of zero also in the definitions of the integrals in (5.2). 

Remark 13. The solution to the linear differential equation (5.3) is the sum over 
the set of solutions with inhomogeneous terms p(x) which are powers of x. These are 
Bernoulli differential equations, so the standard substitution f = q 1 esp reduces 
the problem to the differential equations 

(5.5) + i(l - degp)(f> , (x)f(x) = 1 - degp. 
ax 

In that sense, the only form to be studied is the one with constant left hand side in 
(5.3). This route has not been followed here, because the connection to the initial 
value problem lim^-yo Q — ► is lost by changing the variable to f which is an inverse 
power of q. 

Remark 14. This is an analytic variant of Levin's collocation method of casting 
an approximate solution [35, 13]. It is geared to a power-series solution (the degree 
of the polynomials q increases with their index) and in that sense another method 
to generate the Taylor series at x = 0. 

5.2. Asymptotics. The asymptotic behavior 

(5.6) *~E| = EM^' 

j'=i j>i 

is found by substituting X = x^ 1 in the polynomials p and <j)' , replacing dq/dx = 
—X 2 dq/dX in (5.3), multiplication of the equation by X dc6 ^ and solving the dif- 
ferential equation for q(X) with a power series ansatz [44], see Appendix B. The 
leading term of the series has order X des< ^ _do sP, where we may assume that this 
series exists by enforcing dcg0' > degp as discussed in Remark 4. 

Example 2. If p = x 2 and <f> = x + x A , <j)' = 1 + 4a; 3 , the differential equation is 
X- 2 = iq(l + AX- 3 )-X 2 q'(X), or X = iq(X 3 + 4) - X 5 q'(X), which is solved by 
i i 1 i 
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exact 
10 
20 
30 
40 




-0.4 - w 

-0.5 I ' 1 1 1 — 1 1 

0.5 1 1.5 2 2.5 3 

x 

Figure 6. The exact value of function of upper limit x 

for the example p = x 2 , (f> = x + x A , and four polynomial approxi- 
mations that gather 10, 20, 30 or 40 orders of the expansion (5.2). 
The approximations stay on the solution for approximately 1 to 3 
oscillations and exit abruptly. 



The efficiency of this series to estimate the amplitude of the oscillations of I p ^ for 
large x is illustrated in Figure 7. 

5.3. Anticyclic. So far, the expansion (5.2) has anchored the values of I p ^{x) at 
x = 0, and the asymptotic expansion (5.6) has anchored them at x = oo. The two 
solutions q of different initial conditions miss each other by a gap I p ^(po) 1 

(5.8) e 1 ^ ]T q-ik+i = J P ,*(oo) + e***) £ h 3 XK 

k>0 j>l 

The solutions differ by a multiple of I p ^{oo)e~' 1 ^ , the solution to the homogeneous 
differential equation, following the general theory of linear differential equations. 

Remark 15. I p ^(oo) might be called the "complete'' generalized Fresnel Integral. 

If we would start the perturbative accumulation of q not with the estimate p as 
in (5.2) but with some constant qo = e. the repeated integration and multiplica- 
tion with <f>' in (5.2) would generate the term —i J ecj)'dx = —iecj), then the term 
—i j '(—ie<j))(f>' dx = e(—i<p) 2 /2, and eventually the entire general solution ee~ l< ^ of 
the homogeneous differential equation. The equivalent procedure applied to (5.2), 
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0.5 



0.4 



0.3 



E 0.2 



0.1 



-0.1 



0.5 



1.5 

x 



~i r - 

exact 
order X 
order X* 
order X 




2.5 



Figure 7. The exact value of as a function of upper limit 

x = l/X for the polynomials p and (j) °f Example 2, plus three 
residual asymptotically horizontal curves 3[7 - (~iX / '4)], 



Wp,4 



(-iX/A + iX 4 /16)} and Q[I p ^ - {-iX / A 



iX 4 /16 — X 5 /16)] after up to three leading components of the 
asymptotic have been removed. 



splitting <72fc-i o,t each iteration into a multiple of <fi and remainder, yields a rep- 
resentation q = ^2 S>0 q s (x)(i<p) s / si, where q s are polynomials of degree less than 
deg0. 

The effect on q is that ^<?2fe+i picks up an ant i- resonant, counter-cycling phase 
7p j 0(oo)e -4 ^^^ as x increases, mentioned in Remark 8 and observed again in Figure 
6.' 

The calculation of I(oo) is accomplished by Taylor expansion of all contributions 
of (j>, setting apart the term of highest polynomial degree, using (3.9). 

Algorithm 3. Let (f) = ^2 l j=i a j x ^ > ' = deg0. Obtain Taylor series coefficients t 

by letting p(x)e l ^j =1 ajX — ^2j> tjX J — optionally with an intermediate expansion 
of the exponential factor in complete exponential Bell polynomials [34] . Then 



(5.9) 



j>0 



(1+3)/' 
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Results of such calculation are put into Table 4. Cases that can be reduced with 
Remark 4 or Remark 11 or Eq. (1.14) have been excised. 
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Table 4. Examples of real and imaginary parts of integrals 
at infinite upper limit. 
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1 
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-0 


08125607922853098 





12268865661476232 


1 
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X 
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1 
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X 
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2x 
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-0 
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Unfortunately some fiuky maneouvre is still needed to decide on a transition 
value that divides the x-interval into the low region of good convergence of the sum 
^2P2k+i and the high region of good convergence of the sum ^ hjX J , depending 
on requirements on accuracy and number of terms of both series to be computed. 
In practise one picks some x larger than 1 where c^fc+i converges, computes the 
asymptotic series at the associate X = 1/x, estimates the error in the asymptotic 
part (which adds some artistic flavor to the concept) and obtains a point spot 
estimate of /^(oo) via (5.8) subtracting both and multiplying with e l ^ x \ 

Example 3. Continuing Example 2, one may compute a stable value of p(x) near 
x pa 1.8 with 40 terms, illustrated in Figure 6, and compute an associate p(X) 
again near x ~ 1.8 with the first three terms illustrated by the collocation of the 
X 4 and X 5 curves in Figure 7. The result is ip^oo) ~ —0.07116+ 0.20606z with 
approximately two valid digits compared to Table 4- The imaginary part matches 
the step height as x — >■ oo in Figure 7. 



The integral over the product of a polynomial p(x) by a phasor cxjp[i<p(x)] defined 
by another polynomial (f>(x) has been reduced to Incomplete Gamma Functions if p 
and <p are powers of x, and to power series in x and asymptotic series in 1/x for the 
generic case which differ by the general solution of an associated linear first-order 
differential equation. 



Some elementary examples of series reversions (4.2) are gathered in Table 5. It 
shows the polynomial in x, the initial few coefficients f3j of its reversion, and — if 
available — the label of these coefficients in the OEIS [41]. 



6. Summary 



Appendix A. Reversion of basic polynomials 



Table 5: Basic Examples of Scries reversions 



x + 2x 2 



9 ^ 

X — x — X 



x — x 2 + x 3 



X + X 5 



X + x 2 + X 3 



X + X 3 



y = •• 

X + X 2 



X 



X 



.2 



ft, [41] 

1, -1, 2, -5, 14, -42, 132, -429, 1430, -4862, 16796, A168491 

-58786, 208012, ... 

1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, A000108 
208012, ... 

1, 0, -1, 0, 3, 0, -12, 0, 55, 0, -273, 0, 1428, 0, -7752, 0, A001764 
43263, 0, -246675, ... 

1, 0, 0, -1, 0, 0, 4, 0, 0, -22, 0, 0, 140, 0, 0, -969, 0, 0, A002293 
7084, 0, 0, -53820, ... 

1, -2, 8, -40, 224, -1344, 8448, -54912, 366080, A052701 
-2489344, 17199104, ... 

1, 0, 0, 0, -1, 0, 0, 0, 5, 0, 0, 0, -35, 0, 0, 0, 285, 0, 0, 0, A002294 
-2530, 0, 0, ... 

1, -1, 1, 0, -4, 14, -30, 33, 55, -429, 1365, -2652, 1428, A103779 
12920, -64600, ... 

1, 1, 1, 0, -4, -14, -30, -33, 55, 429, 1365, 2652, 1428, A103779 
-12920, -64600, ... 

1, 1, 3, 10, 38, 154, 654, 2871, 12925, 59345, 276835, A001002 
1308320, ... 
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Table 5: Basic Examples of Scries reversions (cont.) 

» = ••■ 01,02,03,... [41] 

z + x 6 1, 0, 0, 0, 0, -1, 0, 0, 0, 0, 6, 0, 0, 0, 0, -51, 0, 0, 0, 0, A002295 

506, 0, 0, 0, ... 

x + 2x 3 1,0, -2, 0, 12, 0, -96, 0, 880, 0, -8736, 0, 91392, 0, A153231 

-992256, 0, 11075328, ... 
x + x 2 +x 4 1,-1,2,-6,20,-70,256,-969,3762,-14894,59904, A049140 

-244088, 1005452, ... 

x + 3x 2 1, -3, 18, -135, 1134, -10206, 96228, -938223, 9382230, A005159 

-95698746, ... 

x + x 7 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, -70, 0, 0, A002296 

0, 0, 0, 819, ... 

1, 0, -1, -1, 3, 7, -8, -45, 0, 264, 273, -1365, -3192, A217359 
5508, 27132, -7752, ... 

1, -1, 2, -5, 13, -35, 96, -264, 720, -1925, 4966, -12038, A063028 
25907, -41310, ... 

x-x 2 -x 5 1, 1, 2, 5, 15, 49, 168, 594, 2150, 7931, 29718, 112814, A063021 

432957, ... 

x + 2x 2 +x 3 1,-2,7, -30, 143, -728, 3876, -21318, 120175, -690690, A006013 

4032015, -23841480, ... 
x + x 8 1, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 8, 0, 0, 0, 0, 0, 0, A007556 

-92, 0, 0, 0, ... 

x + 2x A 1, 0, 0, -2, 0, 0, 16, 0, 0, -176, 0, 0, 2240, 0, 0, -31008, A217360 

0, 0, 453376, ... 

x + x 3 +x 5 1, 0, -1, 0, 2, 0, -4, 0, 5, 0, 13, 0, -147, 0, 816, 0, -3534, 

0, 12650, 0, ... 

x + x 2 +x 6 1, -1, 2, -5, 14, -43, 140, -474, 1650, -5863, 21170, 

-77441, 286372, ... 

x + x 2 + 2x 3 1,-1,0,5,-16,14,96,-495,880,2002,-17888,48178, A217361 

19040, -665380, ... 

x + 2x 2 +x 4 1, -2, 8, -41, 236, -1456, 9412, -62922, 431464, A217362 

-3017894, 21448024, ... 
x + Ax 2 1, -4, 32, -320, 3584, -43008, 540672, -7028736, A052704 

93716480, -1274544128, ... 
x + x 9 1, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 9, 0, 0, 0, 0, 0, A062994 

0, 0, -117, ... 

x + x i + x 5 1, 0, 0, -1, -1, 0, 4, 9, 5, -22, -78, -91, 105, 680, 1224, 

0, -5700, -14630, ... 
x + x 3 +x 6 1, 0, -1, 0, 3, -1, -12, 9, 55, -66, -267, 455, 1323, 

-3060, -6528, 20298, ... 
x + x 2 +x 7 1, -1, 2, -5, 14, -42, 131, -420, 1375, -4576, 15431, 

-52598, 180887, ... 

x + 3x 3 1, 0, -3, 0, 27, 0, -324, 0, 4455, 0, -66339, 0, 1041012, 0, A217363 

-16953624, ... 

x + x 2 +x 3 + x A 1, -1, 1, -1, 2, -7, 22, -57, 132, -308, 793, -2223, 6328, A063019 

-17578, 47804, ... 

x + 2x 2 +x 5 1, -2, 8, -40, 223, -1330, 8304, -53592, 354645, 

-2393380, 16409224, ... 
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Table 5: Basic Examples of Scries reversions (cont.) 

y = --- [41] 

x + 3x 2 +x 3 1, -3, 17, -120, 948, -8022, 71106, -651717, 6126175, 

-58736535, 572178165, ... 
x + x w 1, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, A062744 

0, 0, 0, ... 

x + 2 x 5 1, 0, 0, 0, -2, 0, 0, 0, 20, 0, 0, 0, -280, 0, 0, 0, 4560, 0, 0, A217364 

0, -80960, ... 

x + x A + x 6 1, 0, 0, -1, 0, -1, 4, 0, 10, -22, 6, -91, 140, -120, 816, 

-1020, 1710, -7315, ... 
x + x 3 + x 7 1, 0, -1, 0, 3, 0, -13, 0, 65, 0, -351, 0, 1995, 0, -11764, 

0, 71307, 0, ... 

x + x 2 +x 8 1, -1, 2, -5, 14, -42, 132, -430, 1440, -4928, 17160, 

-60606, 216580, ... 

x + 2 x 3 + x A 1,0, -2, -1, 12, 14, -92, -180, 770, 2266, -6552, -28210, 

53452, 347616, ... 

x + x 2 + 2 x 4 1, -1, 2, -7, 26, -98, 388, -1599, 6754, -29062, 127036, 

-562614, 2518852, ... 
x + x 2 + x 3 + x 5 1, -1, 1, 0, -5, 21, -58, 108, -50, -638, 3471, -11323, 

24843, -20910, -123216, ... 
x + 2x 2 +x 6 1, -2, 8, -40, 224, -1345, 8464, -55092, 367840, 

-2505360, 17338886, ... 
x + 2 x 2 + 2 x 3 1, -2, 6, -20, 68, -224, 672, -1584, 880, 22880, -215072, 

1414400, -8012032, ... 
x + 3x 2 + x 4 1, -3, 18, -136, 1152, -10458, 99472, -978453, 9871686, 

-101590654, ... 

x + 5x 2 1, -5, 50, -625, 8750, -131250, 2062500, -33515625, A156058 

558593750, -9496093750, ... 
x + x 11 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 0, 

0, 0, 0, ... 

x + x 5 + x 6 1, 0, 0, 0, -1, -1, 0, 0, 5, 11, 6, 0, -35, -120, -136, -51, 

285, 1330, 2310, ... 

x + x 4 + x 7 1, 0, 0, -1, 0, 0, 3, 0, 0, -11, 0, 0, 42, 0, 0, -153, 0, 0, 

469, 0, 0, -690, ... 

x + x 3 + x 8 1, 0, -1, 0, 3, 0, -12, -1, 55, 11, -273, -91, 1428, 680, 

-7744, -4845, 43092, ... 
x + x 2 +x 9 1, -1, 2, -5, 14, -42, 132, -429, 1429, -4851, 16718, 

-58331, 205632, ... 

x + x 3 + 2 x A 1, 0, -1, -2, 3, 14, 4, -90, -165, 396, 1911, 0, -15372, 

-25704, 85272, ... 

x + 2x 3 +x b 1, 0, -2, 0, 11, 0, -80, 0, 665, 0, -5980, 0, 56637, 0, A118969 

-556512, 0, 5620485, ... 

1, -1, 2, -6, 19, -63, 220, -795, 2942, -11099, 42536, A063030 
-165126, 647955, ... 

1, 1, 3, 10, 39, 161, 698, 3126, 14360, 67276, 320229, A063022 
1544257, ... 

1, 1, 2, 6, 21, 77, 292, 1143, 4592, 18821, 78364, 330512, A063023 
1409149, ... 
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Table 5: Basic Examples of Scries reversions (cont.) 

» = ••■ 01,02,03,... [41] 

x + x 2 + x 3 + x 6 1, -1, 1, 0, -4, 13, -22, -3, 165, -638, 1371, -832, 

-6949, 34850, -89216, ... 
x + 2x 2 +x 7 1, -2, 8, -40, 224, -1344, 8447, -54894, 365860, 

-2487056, 17177264, ... 
x + x 2 +3x 3 1, -1, -1, 10, -22, -42, 438, -1023, -2695, 26455, 

-60073, -203320, 1845928, ... 
x + 2x 2 +x 3 +x 4 1, -2, 7, -31, 155, -833, 4696, -27393, 163944, 

-1001022, 6211049, ... 
x + 3 x 2 + x 5 1, -3, 18, -135, 1133, -10185, 95904, -933768, 9324320, 

-94969215, ... 

x + Ax 2 +x 3 1, -4, 31, -300, 3251, -37744, 459060, -5773548, 

74474455, -979872036, ... 
x + 6x 2 1, -6, 72, -1080, 18144, -326592, 6158592, -120092544, A156128 

2401850880, ... 

x + x 2 + x 3 + x 4 + x 5 1, -1, 1, -1, 1, 0, -6, 27, -83, 209, -455, 845, -1169, A217365 
272, 5916, -29070, ... 

To support economical evaluation of the reversions, the subsequent list shows 
(homogeneous) hypergeometric recurrences in annihilator notation. 

These formats occur because the /? are determinantal forms of the a-cocfficients 
of the polynomials [19, 28, 31]. From the point of view of Lagrange's inversion 
formula, 

At-l ( „ \ ' At-\ | 

(A.l) t\Pt 

u 

accompanying (4.1), evaluated with Faa di Bruno's formula [18, 15, 32, 45], it 
portraits that the number of nonzero higher order derivatives of the inner function 
(the polynomial in the denominator) is limited by the order I. Expansion of the 
determinant along the (sparse) final row of di-Bruno's determinant establishes the 
recurrence. 

Two-term hypergeometric recurrences in this list lead to representations of (3 as 
product-ratios of T- values, and to generating functions for the /3-sequence of the 
hypergeometric form. 

The simplest group are the cases y — x + kx* with integer k and t where the 
coefficients (3 are generalized Catalan numbers [14]: 

• The coefficients for series reversion of y = x + kx 2 are essentially products 
of Catalan numbers with powers of k and obey 2-term recurrences. Solving 
the quadratic equation for a; yields [9] 

(A.2) x = [(1 + Akyf' 2 - l]/(2fc) = Y, (H 2 ) (4M7(2fc) 



XLi a ^ xk J \ x=0 dxt 1 (EUk) *k+ix k y , B=0 



!>1 



£f 2 /:, 2 ^'-'<-»)<- 



i>i v 

So (3 3 = W-fly-ky- 1 and jfy + 2k{2j - 3)/?,_i = 0. 
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• The non- vanishing coefficients of y = x + x 3 obey (2k — l)(2fc — 2) 7 fc + 
3(3k - 4)(3fc - 5) 7 fe-i = 0, where lk = fek-i- 

• The non- vanishing coefficients of y = x + x A obey (3fc — 3)(3/c — 4)(3fc — 

2) 7fc + 4(4fc - 5) (4* - 7)(4fc - 6) 7fe -i = 0, where 7fe = 3fc _ 2 . 

• The non- vanishing coefficients of y = x + x 5 obey (4fc — 4)(4fc — 5)(4fc — 

3) (4fc-6) 7fe + 5(5fc-9)(5fe-8)(5fe-7)(5fc-6) 7fe _i = 0, where lk = /3 4 fe-3- 

• The non-vanishing coefficients of y = x + 2x 3 obey (2k — l)(k — l)jk + 
3(3fc - 4)(3fc - 5) 7fc -i = 0, where 7fe = /? 2 fc-i- 

• The non- vanishing coefficients of y = x + 2x A obey (3fc — 3)(3fc — 4)(3/c — 
2) 7fe + 8(4fc - 5)(4fc - 7) (4* - 6) 7fe _i = 0, where 7fc = /? 3fc - 2 . 

Other recurrences are 

• The coefficients of y = x + x 2 + x 4 obey 1147j(j - l)(j - 2)(5 3 + 8(j - 

1) (j - 2)(647j - 738)/3j_i + 4(j - 2)(224j 2 + 1504 j - 5157)/%_ 2 + 8(800j 3 - 
5040j 2 + 8746j - 2655)ft_ 3 - 192(4j - 15)(2j - 7)(4j - 17)/3j-_ 4 = 0. 

• The coefficients of y = x + x 3 + x 4 obey 124j(j - l)(j - 2)f3j + (j - l)(j - 

2) (7j - 88)/3j_i + (j - 2)(870j 2 - 3456j + 3347)/3,-_ 2 + (1243j 3 - 9870j 2 + 
25869j - 22490)^-3 + 8(4j - 15)(2j - 7)(4j - 17)^_ 4 = 0. 

• The coefficients of y = x + 2x 2 + x 3 = x(l + x) 2 obey 2j(j - l)/3j + 3(3j - 

2) (3j - 4)&_i = 0, so = IlU(-3)(3i - 2)(3i - 4)/[2i(* - 1)]. 

• The coefficients of y = x + x 2 + 2x 3 obey 7j(j - l)/3j + 16(j - l)(2j - 

3) ^-1 + 12(3j - l)(3j - 7)&_ 2 = 0. 

Appendix B. Laurent Series 

The coefficients of the Laurent expansion (5.6) are obtained by inserting the 
series and the polynomial <j> (1.13) into (5.3): 

i 

(B.l) X~ m = lq (X)Y J J^ 3 X^ - X 2 q'(X), 

3 = 1 

where we assume that only the forms p = x m = 1/1'", m < I — 1, need to be 
studied because the differential equation is linear and because higher powers m 
have been decimated as proposed in Remark 4. We compare coefficients of equal 
powers of X in 

oo Z oo 

(B.2) X~ m = i^2h t X t ^2ja j X~ j - X 2 ^2th t X t -\ 

t=i j=i t=i 

Elementary resummation of the double sum on the right hand side suggests 

/ oo I oo 

(B.3) X~ m = iJ2 J2 i'^'-i^'"^^-- 1 ' Y. is '.^ i- v ~- 

S=— I j= — 3+1 S—lj — 1 S—l 

For s < 0, only the first double sum on the right hand side contributes, so the ht 
are obtained recursively for increasing t = 1, 2, . . . , I via 

-l 

(B.4) ilaiht = 6 t -i,- m - i ^ (j + l)a J+ ih t+j , l<t<l. 

j=-t+i 

For s > 0, the left hand side of (B.3) vanishes and only the second and third sum on 
the right hand side contribute, so the ht are obtained recursively for t = l+l, 1+2, . . . 
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via 

-l 

(B.5) ilaiht = (t - I - l)h t -i-i - i (j + l)a j+ ih t +j, I < t. 

j=-l+i 

with the convention that ht = whenever t < 0. 
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